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ABSTRACT 

Aims. We derive lists of proper-motions and kinematic membership probabilities for 49 open clusters and possible open clusters in 
the zone of the Bordeaux PM2000 proper motion catalogue (-1-11° < 6 < +18°). We test different parametrisations of the proper 
motion and position distribution functions and select the most successful one. In the light of those results, we analyse some objects 
individually. 

Methods. The segregation between cluster and field member stars, and the assignment of membership probabilities, is accomplished 
by applying a new and fully automated method based on both parametrisations of the proper motion and position distribution func- 
tions, and genetic algorithm optimization heuristics associated with a derivative-based hill climbing algorithm for the likelihood 
optimization. 

Results. We present a catalogue comprising kinematic parameters and associated membership probability lists for 49 open clusters 
and possible open clusters in the Bordeaux PM2000 catalogue region. We note that this is the first determination of proper motions for 
five open clusters. We confirm the non-existence of two kinematic populations in the region of 15 previously suspected non-existent 
objects. 

Key words. Galaxy: open clusters and associations. Methods: data analysis. Methods: statistical 



1. Introduction 

Once the problem of selecting their physical members is re- 
solved, open clusters are widely respected to be a most valu- 
able tool for undertaking studies of our Galaxy and stellar as- 
trophysics. These objects have been used, for example, to de- 
termine the spiral structure of the Galaxy and investigate star 
formation and evolution process es. They are particularly impor- 
tant as tracers of the dynamics d Frinchabov & Mai ewskil 12008^ 
and the chemical evolution of our Galaxy's disk (Friel 1995). In 
the advent of high pre cision astrophysical surveys such as Gaia 
dPerrvman et al.ll2001h . their contribution to astrophysical stud- 
ies should become increasingly important. 

However, to establish a coherent understanding of our 
Galaxy, one needs to use a significant number of open clusters 
with consistently measured astrophysical para meters, and there 
have been numerous efforts in this direction (Kharchenko et al.' 
120051; iBragaglia & Tbsil 120061: iF rinchabov & Maiewski 2008). 
However, as noted by'Frinchabov & Maiewski (2008), the inac- 
curacy of cuiTent membership determinations poses difficulties 
in conducting these studies on a large scale. 

The main advantage of using open clusters in these studies 
is that once the complex problem of stellar membership is re- 
solved one can derive their main physical parameters such as 
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distance, age, and metallicity. This membership determination is 
traditionally performed using the stellar kinematics, but in prin- 
ciple one could use a multidimensional space, using for exam- 
ple, spatial and kinematic information, as in Z hao et alJ 
or CMD-isochrone information, as in KJiarchenko et al.l 
Nonetheless, we argue that, when analysing these objects, one 
should rely more on kinematics and as little as possible on a 
single age CMD-based model analysis as in the aforementioned 
study. This is because, as it has been reported in the litera- 
ture, the star formation in some open clusters could be non- 
coeval, as seems to be true in NGC 3603 dEisenhauer et alJl998h 
and 14 other open clusters dStrobe]|ll992ir An indication of an 
abundance sp read has been reported for one of these objects 
dFrinchabov e t al. 2008). 

Based on these results, we undertook a purely kinematic de- 
termination of open cluster membership probabilities for objects 
located in the zone covered by the high precision proper mo- 
tion catalogue PM2000 ( Ducourant et al 2006) . We used a fully 

automated optimization method and a set of modified parametri- 
sations for the probability distribution functions based on Zhao 
et al. (1990, 2006). 

This paper is organised as follows. In Sect. 2, we describe the 
data we used and its selection process. In Sect. 3, we present the 
methodology and algorithms chosen to obtain the membership 
lists and the cluster kinematic parameters. In Sect. 4, we describe 
the validation of the method. In Sect. 5 we present our results, 
and our comment on some individual objects. Finally, in Sect. 6 
we present the conclusions of our current study. 
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2. Data 

2.1. The PM2000 catalogue 



The PM2000 catalogue dPucourant et alj|2006l) is a proper mo- 
tion catalogue that comprises about 2.6 million stars in the dec- 
lination zone +11° < 6 < +18° and contains positions and 
proper motions on the ICRS (International Celestial Reference 
System), as well as meridian magnitudes Vm- It was derived 
from the compilation of systematic drift-scan observations in 
the Bordeaux Carte du Ciel Zone with th e Bordeaux auto- 
mated meridian circle dViateau et al.l 1 19991) carried out over 
four years, the reduction of 51 2 Carte du Ciel plate s (epoch 
t X 1900) of the Bordeaux zone jRapaport et alj|200lT) scanned 
at the APM Cambridge, and the catalogues AC2000.2 (f « 
1907), USNO-A2.0 (f ~ 1950) and the unpubhshed USNO 
Yellow Sky (YS3, t x 1978). The positional precision ranges 
from 50 to 70 mas, while the proper-motion precision varies 
from 1.5 to 6 mas yr depending on the magnitude. All the 
data was analysed using a global iterative astrometric reduc- 
tion (Tei xeira et al. 1992; ' Benevides-Soares & Teixeiral Il992t 
iDucourant & Rapaport 199iIk 

The catalogue is complete to Vm - 15.4, with a limit- 
ing magnitude of Vm = 16.2, and typical error of 0.03 mag 
(9.5 < Vm < 13.5). In addition, a cross identification bet ween 
all sources in the PM2000 and 2MASS (ICutri et al.ll2003h was 
performed, so the PM2000 catalogue also includes 2MASS pho- 
tometry information for its objects. 

2.2. The clusters sample 

The starting point of our analysis is a list of 49 open clusters 
inside the PM200 declination zone found in the D07 catalogue 
jDias et alj2002al) . We performed a vis ual inspection of all these 
clusters by using the Aladin Sky Atlas ('B onnarel et al.ll2b00i) to 
verify their coordinates. During this visual check, we noticed 
that the clusters Berkeley 29, 43, 45, and 47 needed to be slightly 
recentered about 3'25", 3'39", ri3", and r27", respectively. 

The objects were then separated into four different classes: 
A-reference, B-known, C-known without proper-motion determi- 
nation, and D-others (including doubtful objects). An object was 
classified as D and not C when one or more of the following con- 
ditions applied: 

1 . No entry in the WEBDA database; 

2. Classification as not found, dubious, no cluster, or non- 
existent NGC in the D07 catalogue; 

3. Existence of some study that excludes the physical clustering 
of its stars (such as for NGC 1807). 

For each cluster, we extracted from the PM2000 catalogue 
all the stars inside a circular area around the apparent cluster 
centre. We multiplied the cluster diameter by 1 .5, to include most 
of the cluster members in the extracted zones. Two arcminutes 
were then added to ensure that even for a small cluster, located 
in a poorly populated region of the Galaxy, a sufficient number 
of field stars would be available in the extracted zone to enable 
us to discriminate between the cluster and the field (as in the 
case of NGC 7772). We visually inspected the DSS images of 
all the extracted zones, to ensure that the objects were inside the 
extracted regions if they were clearly identified. We note that this 
extended region could add some noise (in the form of field stars) 
to a small object located in a densely populated region. All the 
extractions were performed by using automated PERL scripts 



and CDS's vizquery tool. Table[T]shows the input Ust of clusters 
that we used, as well as the class that we assigned to each object. 

The data used to obtain the cluster kinematical parameters 
were selected from the extracted data by rejecting on the basis 
of proper motion errors (e^ < 6.0 mas yr ') and the total proper 
motion (|yu| < 30.0 mas yr"'). This was justified since high /i 
field stars cause the proper motion distribution to be flattened, as 
previously noted in the literature (Bal aguer- Nunez et al, 2004bl 
and references therein). Nonetheless, we computed membership 
probabilities for all the stars in the extracted zones. 

3. Methods 

3.1. Mathematical model 

The traditional way of conducting membership assignment orig- 
inates in the seminal works of IVasilevskis et al.l (11958 ^ and 
Sanders (1971). Those two works provided the basic ideas for 
developing parametric membership analysis. The derived math- 
ematical model is based on the assumptions that there are two 
distinct kinematic populations in the observational field of the 
cluster and that the proper motion distributions of those two pop- 
ulations can be parametrised by bivariate Gaussians. 

In these studies, an elliptical function was used to describe 
the field population's proper motion, while a circular one was 
used for the cluster. Nonetheless, if either an external gravita- 
tional influence or tidal effects were to act on different parts of 
the cluster, the cluster's proper motion distribution could be sig- 
nificantly affected, causing it to deviate from a circular function. 
For this reason, we tested four parametrisations of the probabil- 
ity distribution functions (PDF) in this study. 

We adopted three variations of the general form of the PDF 
established by Zhao & He ( 1990) to take account of the obser- 
vational errors in each individual point. The proper-motion dis- 
persion parameters obtained from these PDFs are therefore the 
cluster and field intrinsic ones, which are independent of the ob- 
servational errors in the proper motions. The first variation is a 
circular distribution, which is the most accurate representation 
of the PDF for a non-disturbed object with symmetric observa- 
tion errors. The second variation is elliptical (allowing for the 
correlation coefficient) and the last a PDF in which we assume 
that the intrinsic dispersion can not be observed because of th e 
size of the errors (as adopted by iBalaguer-Nunez et"al ] l2004bh . 
We also tested a fourth PDF in which the {a , S) positions of th e 
individual stars are taken into account, as in IZhao et alJ (l2006l) . 
Nonetheless, unlike this study, which used the radial distance of 
the star to the centre of the cluster, we directly used the stars' co- 
ordinates and considered the cluster centre as a free parameter. 

The adopted PDFs can therefore be written as a mixture of 
proper-motion (<1>) and position (*!*) PDFs 

^^(bc'V. + fbf^f, (1) 

where c and / correspond to cluster and field parameters. 

The proper motion PDFs, Oc and <l>/, are assumed to be 
Gaussians of the form 
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While the and *P/ functions depend on whether we take into 
account the position of the stars or not. We refer to as 2D imple- 
mentations those that do not take into account the positions, and 
4D those that do. These parametrisations are given by 



1 



1+^ ' exp 



{ 2(1 
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In Eq.|2] above, fiaj, fJ-sj represents the proper motion of the 
;-th star, e^^ , e^,, the proper motion errors of the i-th star,/v„, jjg 
the mean proper motions of the cluster and field stars (depending 
on the index of O), cr^, , cr^_, the intrinsic proper-motion disper- 
sions of the cluster and the field, and p represents the correlation 
coefficients of the cluster and field proper motion distributions. 
In Eqs.[3]and|4] rtc, 1 - Wc are the fraction of the cluster and field 
stars, g is the ratio of the cluster to the field stars, a,-, 5, are the 
positions of the i-th star, a, 6 are the position of the cluster, cr„, 
(Ts represent the dispersion, and ppos the correlation coefficient 
of the position distribution. 

The parametrisations presented above leave us with a vec- 
tor of 11 or 16 parameters to be determined, depending on 
whether we consider 2D or 4D distributions. 



3.2. Model parameter estimation 

During the construction of the PM2000 catalogue, a global 
method was used to solve the star's astrometric parameters by 
selecting mean epochs such that the catalogue covariance matri- 
ces are diagonal, and the parameters obtained can be considered 
independent and identically distributed (iid). We are therefore 
allowed to write the likelihood function for a given region with 
N stars as 



(5) 



Following the maximum likelihood principle the most probable 
is that for which X takes its maximum value. A monotone 
transformation does not affect a maximum, and we can therefore 
write - arg © max X. - arg © max X, where 



■£ = ^ log O, (a/, 6i,Haj,ii6j, 0) 



(6) 



Usually the is found to be the solution of a non-linear system 
of equations constructed from dX./d & = 0, in which some sort 
of iterative proce dure is used (as in ISanderslll97lt IZhao & Hj 
[T990tlDiaset"al ] |2OO60 . However, this approach depends on the 
initial chosen value of each unknown parameter. Therefore, the 
solution to the problem will converge to the real, physical one, 
provided that there are no local optima, or the initial value is not 
far from the physical optimum. We note that the last condition 
may be cluster dependent, and for many objects, we will indeed 
have local optima. Thus, one cannot guarantee that the obtained 
solution corresponds to the physical one. 

To help us obtain global optima, we chose to use a strat- 
egy based on evolutionary computing to solve the optimization 
problem, the genetic algorithm (GA). This is an adaptive search 



heuristic for finding optimal solutions, which has been used 
several times in astronomy (iHetem & Gregorio-H etemI 120071 : 
Howley et al. 2008). However, the number of generations neces- 
sary to attain convergence can be quite large, since GAs are gen- 
erally slow to converge. On the other hand, these algorithms are 
usually very good at finding the region of global maximum. In 
this study, we decided to use an optimiser available in the R envi- 
ronment (R Development Core Team 2008) called RGENOUD, 
which is a mixture of a GA and a derivative-based hill-climbing 
algorithm. The GA identifies the region of the global maxi- 
mum, while the quasi-Newton hill-climbing algorithm Broyden- 
Fletcher-Goldfarb-Shanno (or BFGS) based on the fitness func- 
tion's derivatives, is used to optimise the best (higher likelihood) 
solutions and, therefore, find the region's maximum (which may 
be a global maximum). For a complete description, we refer to 
Mebane & Sekhon (2007). 

In principle we could allow the algorithm to search for a so- 
lution in R" or R'^ freely. Nonetheless, we know that, physi- 
cally, the search region can be constrained for some of the pa- 
rameters, such that if an optimal point were to exist outside this 
region, it would not have any physical meaning. For example, 
if most cluster stars were inside the extracted region, we could 
assume that the position of the cluster center would never be 
smaller or greater in value than the positions of the most extreme 
stars in that region. The minimum and maximum of the cluster 
and field proper-motion components should never be smaller or 
greater than the minimum and maximum proper-motion com- 
ponents of the stars themselves. We should also ensure that the 
correlation coefficient moduli never equals unity, otherwise we 
could encounter problems related to divisions by zero. We also 
constrain the cluster's proper motion dispersions to have a max- 
imum of 5 mas yr"', which is an estimated upper limit to the 
average errors in the individual proper motions. Finally, we al- 
low the field's proper motion (and cluster's position) dispersions 
to be almost free, adopting a very high limit of six times the 
standard deviation of the individual proper motions (and stellar 
positions). 

Hence, given the four R' sets Mq, ^^^^^^^ of all the stellar posi- 
tions and proper motions in the region, we adopted the following 
constraining conditions for the unknown parameters: 



(7) 



if 2D 
if 4D 



min^M^J < fia,c/f < max(M^„) 
min (m^,) < fis,c/f < max (m^,) 
0.01 < o-^^^ < 5.00 
0.01 <c7-^„, <6-o-(m^„) 
0.01 < o-^,'^ < 5.00 
0.01<C7-^,', <6-tr(M^,) 
0is such that <! -0.99 < p.;, < 0.99 

< He < 1 

0.01 <g< length (M„) 
min (M„) <a< max (M„) 
min (Ms) <6< max (M^) 
0.01 <o-„ <6-cr(MJ 
0.01 <o-a<6-cr(M5) 
-0.99 <ppo,< 0.99 

where a-(Ma^s,fj„,iis) represents the standard deviation of the re- 
spective set of parameters. 

When a solution is found, the Hessian matrix CH) is com- 
puted for the solution parameters. We use this matrix to esti- 
mate the individual parameter's errors in the optimization pro- 
cedure; this is necessary because for a regular problem (in the 
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sense of lWaldlfT949l) . as considered here, when one maximises 
the log-Hkehhood function, the covariance matrix on the param- 
eters is closely approximated by the negativ e of the i nver se of the 
Hessian matrix at the solution point (see iPawitanI 120011) . Thus, 
the vector of the variances in the parameters is estimated to be 

0,„,^diag (-•?/-'), (8) 

and we use the square root of the elements of this vector to esti- 
mate the fitting procedure errors. 

Once the cluster and field PDF's unknowns are found, the 
membership problem is solved because we are then able to com- 
pute Pi, the cluster membership probability of the /-th star 

The number of apparent cluster members can be estimated from 
the mixture proportion Uc in the following way. First, from the 
list of stars used during the optimization process (the list with the 
quality cuts explained in Sect. 2.2), a sublist of the Nopt*nc most 
probable stars (where Nopt) is created. From this list of most 
probable stars, the probability of the least probable star P„,„, is 
obtained. Finally, the number of members in the extracted re- 
gion (without any probability cut applied) can be estimated from 
the number of stars in that region with membership probabilities 
greater than Pmin- 



4. Validation 

We verify the reliability of our work in four steps. We firstly vali- 
date the automatic optimization method using very precise data. 
Secondly, we test the new 4D parametrisation by comparing it 
with another similar parametrisation and method used before in 
the literature. Thirdly, we validate the data by comparing our 
results for a well known cluster with the literature. Finally, we 
perform a comparison of individual members for one "A" and 
one "B" class cluster. 

4.1. Optimisation method 

Very precise proper motion data are availa ble for the open cluster 
NGC1817. These data were presented bv lBalaguer-Nunez et alj 
(l2004bl) (hereafter BN04), in which proper motions were deter- 
mined from 25 plates covering a total time-span of 8 1 years. The 
mean error in the data of more than 80% of the stars is = 1 .55 
mas yr"', while it is 0.97 mas yr"' for 32% of the stars, for ob- 
jects as faint as V < 16. 

Since we use a new optimization method in this work, we 
tested it using these proper motion data. Since we consider the 
method itself, we used the same parametrisation as in BN04, 
and set the intrinsic dispersion of the cluster Gaussian PDF to 
zero. We also applied the same cuts in the data (|yu| < 30 mas 
yr"'). As can be noticed in Table |2] the agreement between our 
results is quite remarkable, indicating that the method adopted 
herein allows us to correctly determine the parameters in a fully 
automatic manner. 

We also performed extensive tests using Monte Carlo ran- 
domly generated sets of Gaussian mixtures. In virtually all cases, 
the parameters of the distributions were recovered correctly. 



Table 2. Method vaHdation forBN04's NGC 1817 data (units of 
H and cr are mas yr"'). 



Work 




fi„ cos 6 


1^6 


^fin COS 6 




P 


Our 


0.252 


0.25 


-0.94 








method 


+0.02 


±0.11 


±0.07 












2.33 


-4.08 


5.11 


5.93 


-0.04 






±0.24 


±0.27 


±0.18 


±0.20 


±0.04 


BN04 


0.261 


0.29 


-0.96 










±0.02 


±0.10 


±0.07 












2.29 


-4.25 


5.69 


6.38 


-0.08 






±0.02 


±0.27 


±0.02 


±0.14 


±0.03 



4.2. 4D parametrisation 

To test the 4D parametrisation, we compared our results for a set 
of open clusters with those obtained using another parametri- 
sation and meth od. The Stochastic Expe ctation Maximization 
(SEM) method dCeleux & Dieboltl Il986l) can be used to de- 
blend the two components corresponding to the field and clus- 
ter, assuming that the data can be described by a mixture of 
4D Gaussians. The SEM algorithm is non-informative and itera- 
tively solves the maximum likelihood equations, with a stochas- 
tic step, for a multivariate mixture of Gaussian distributions. It 
has been adopted to deconvolve t he thin and thick disk popula- 
tions from velocity distributions (Soubiran 1993: Soubir an et aH 
2003). The algorithm initiates at a position decided at random, 
assuming that the sample is a mixture of 2 discrete components. 
We repeat the process 150 times. Most of the time, the algo- 
rithm converges to the same solution. When multiple solutions 
are found, the most frequent one is adopted. 

We used SEM to analyse all the regions on our list. 
Afterward, we compared the solutions with those obtained by 
the 4D parametrisation described herein. The average compati- 
bility between the results, presented in Fig.[T]is very good, with 
a clear distribution centred on (0.25+0.9; -0.46+2.17) mas yr"' 
before a simple 3cr rejection to eliminate outliers, and centred on 
(0.12 + 0.64; -0.13 + 1.45) mas yr"' after the rejection. 

4.3. Data 

To test the data used throughout this work, we applied the 
method described above to the 4D parametrisation to obtain 
the kinematic parameters for the open cluster NGC 2682. This 
object, also called M67, is an old, ~4 Gyr and well known 
open cluster that h as been extensively studied in the literature 
(lYadav et al.ll2008l and references therein). 

By comparing the results obtained for this open cluster with 
SEM and the new method described above, we find that they are 
in close agreement, generally less than the fitting errors com- 
puted for those parameters: A(fia cos 6, fi^) - (0.06,0.02) mas 
yr-i, A(cr^„,os<5,cr^,) = (0.20,0.10) mas yr"', and Ap = 0.02. 

The comparison of the kinematic parameters derived from 
our method with those from previous studies in the literature 
is shown in Table [3] We observe very good agreement between 
all the results within the estimated errors, although we note that 
some of the errors in that table are estimated from the standard 
deviations of member stars, and not from the fitting procedures. 

4.4. Individual membership 

We verify that the members determined by our method are com- 
patible with those found in previous studies by comparing the 
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Table 3. Kinematic parameters determined for NGC 2682. 



CO 
E 

<1 




A|j,aCos6 (mas/yr) 

Fig. 1. Comparison between the result of the fit using the method 
described in this work and SEM. The three filled triangles indi- 
cate points with error bars greater than the plot axis. 



members obtained for our "A" class cluster (NGC 2682) and 
one "B" class cluster (NGC 7036) with previously published 
membership lists. In performing this analysis, we adopted the 
parametrisation that does not take into account the cluster's in- 
ternal dispersion. 

For NGC 2682, 'Yadav et al.J (l2008h (hereafter Y08) com- 
puted the membership probabilities of 2410 stars in its vicin- 
ity. Their study was relatively deep in magnitude, to V~20, al- 
though probably does not probe the entire cluster because of its 
large angular extension. Unfortunately, the authors did not indi- 
cate the number of members found, but a cut at 60% in member- 
ship probability was adopted in their work. For this value, 595 
stars would be considered as members. When analysing PM2000 
data, we obtained 502 members from a total of 1386 stars. We 
note that our catalog is much brighter than that of Y08, but cov- 
ers a more extended area. 

Between those two membership lists, there are 271 stars in 
common. Nonetheless, to perform comparison, we need to take 
into account that only 553 stars are common between our initial 
catalogues, and that among these common objects, 313 are con- 
sidered to be members by ourselves, while 402 are Y08 mem- 
bers. This means that ~87% of our members that could be listed 
as members in Y08 were listed as such, while only ~68% of Y08 
members were listed by ourselves. 

Since PM2000 proper motions are precise in this particu- 
lar field, with 07, < 2 mas yr ' for about 50% of them, but 
in the same magnitude range (V<16) Y08 has only about 28% 
of its stars with such small errors, we attribute this lower value 
of ~68% to a possible contamination of field stars among Y08 
members. 

For NGC 7036. lDias et"al] (|2"006) (hereafter D06) analysed 
69 stars in its vicinity, obtaining 20 members for this possible 
open cluster remnant. We analysed 91 objects close to NGC 
7036, also obtaining 20 members. Adopting the same procedure 
as applied to NGC 2682, we found that ~65% of our members 



Reference 



yu„ cos 6 
(mas yr"') 



(mas yr"') 



This work 

Frinchabov & Maiewski (2008) 
_^as et al. f2002a^ 
Kharchenko et al. (2005') 



-8.32 ± 0.07 
-7.87 ±0.61 
-8.62 ±0.28 
-8.31 ±0.26 



-5.65 ± 0.07 
-5.60 ±0.59 
-6.00 ± 0.28 
-4.81 ±0.22 



that could be listed as members in D06 were listed as such, while 
~68% of D06 members were among our members. 

Taking the results for these two clusters in to consideration, 
we conclude that there is relatively good agreement between the 
members obtained by our method and data with those previously 
published. 



5. Results 

We applied the optimization method and all variations in the 
parametrisation described in Sect. 3.1 to the extracted PM2000 
data. We thus obtained four solutions for each cluster, and have 
to determine the most reliable. We first verify whether each solu- 
tion is in agreement with all others. Next, we visually inspected 
probability histograms, the vector point diagrams (VPD) and star 
charts of members and non-members stars, and the DSS images 
of all clusters, and classified the results as "good", "intermedi- 
ate", or "poor". To avoid any prejudice, we performed a blind 
classification: the names of the objects were not written on the 
diagrams and images while they were being analysed. During 
this check, we noticed that the probability histogram exhibited 
the expected field-cluster distribution with a dual peak distribu- 
tion in most cases. Furthermore, in all cases, the VPD diagrams 
were indicative of accurate fits. Hence, we decided not to use 
those two diagrams when assessing the physical reliabilities of 
the four solutions 

Although this classification of the solutions was subjective, 
we adopted two clear criteria: some kind of structure or central 
concentration of stars in the member star chart (in contrast to an 
almost homogeneous distribution for non-members) and some 
correlation between a certain magnitude range and the member 
stars. The solutions were classified as "good" when a clear con- 
centration of member stars was present (even if there was some 
contamination), and as "intermediate" when the members were 
not so clear concentrated (or when the concentration was very 
off-centred), but there was some correlation between the mem- 
bers and a certain magnitude range. Finally, they were classified 
as "poor" when the members were sparsely distributed without 
any apparent correlation with any magnitude range. 

For NGC 2682, a reference open cluster for which we ob- 
tained "good" results, it is clear from the member diagrams 
that the vast majority of its stars were correctly classified, even 
if the inspection of its non-member diagrams indicates that 
some of this cluster's members were probably not classified as 
such. Moreover, we were able to obtain good results for the 4D 
parametrisation in this cluster, because it is a well-populated 
cluster in the magnitude range of PM2000, which allowed a cor- 
rect determination of its centre and physical dispersion, and its 
kinematic parameters. 

However, of the four parametrisations that we tested, the 
PDF which considers the internal dispersion of the clusters to 
be too small to be observed, provided the least number of results 
classified as "poor". This can be explained physically by the data 
of the distant clusters being affected by this dispersion in the 
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proper motion far less significa ntly than the errors in the indiv id- 
ual measurements, as noted bv Bala guer-Ntinez et alJ (l2004bh . 

After iGunn et aH ([l988), the theoretical prediction for the 
intrinsic velocity dispersion for a Hyades like open cluster in 
a state of dynamical equilibrium is 0.23 km s and the ob- 
served value ba sed on Hipparcos d ata for the Hyades cluster 
is 0.3 km s~' (|de Bruiine et alJ 12001. and references therein). 
iMamaiekl (12010) " reported an intrinsic dispersion of 1 km/s for 
several nearby ob jects, and an up per hmit of 1.1 km s ' in o- 
Persei was set bv iMakarovl ilOOd ). Using the latest proper mo- 
tion catalogues, it would be unrealistic to probe values similar to 
those for any relaxed open cluster that is neither located at very 
nearby distances nor is highly populated: at only a few hundred 
parsecs, the dispersion is already at fias yr ' scales. Thus, the 
use of this additional degree of freedom in the analysis could al- 
low the likelihood function to reach a better fitness value when 
modelling well the field distribution than when segregating field 
from cluster stars. 

We present the results obtained by the PDF that considers 
the internal dispersion of the clusters to be too small to be ob- 
served, for our entire input list in Table|4] It includes all the fitted 
parameters from the cluster and the field PDFs, as well as their 
associated errors. In this table, we also present an estimate of the 
number of apparent members, and the reference class from Table 
1 . This table is divided into three sections, corresponding to the 
quality of the solutions ("good", "intermediate", "poor"). 

In Table |5] (fully available only through CDS), we present 
the membership probability information (computed from Eq. |9]l 
for all the stars in the extracted regions. In addition, we include 
some columns from the PM2000 catalogue, as well as the ob- 
ject's 2MASS identifier when available. 

For some clusters on our list, D 06 also obtained kine matic 
parameters by using UCAC2 data (IZacharias et al.ll2004t) . For 
those objects, we compare our proper motions with those of D06 
work in Fig. |2] We note that the objects tend to have systemat- 
ically positive values of AyU„ cos 6 and Afig alike, which means 
that their proper motions determined by PM2000 data are some- 
what higher than those determined by UCAC2 data, by (0.8, 0.1) 
mas yr The same was found for all the other three parametri- 
sations. 

These results were surprising to us, so we performed several 
tests using the comparison data. During those tests, we found 
that the proper motion differences AyU„ cos 6 were heavily cor- 
related with the cluster's right ascension. This systematic effect 
can be clearly seen in Fig. [3] in which we present the differences 
between the results of the parametrisation that does not take into 
account the cluster's internal dispersion and the D06 catalogue. 

We concluded that the origin of this systematic effect was the 
data itself, since the (a, 6) position of the object was not used in 
3/4 of our parametrisations, it was not used in D06 reduction, and 
the results of all the four parametrisations exhibited exactly the 
same trends. We therefore directly compared the proper motions 
of the common stars between PM2000 and UCAC2 (the proper 
motion catalogue used by D06) inside the cluster regions, and 
noted a clear dependence on the right ascension, which can be 
seen in Fig.lD 

We then cor npared the entire P M2000, UCAC2, U CAC3 
(IZacharias et all 1201 Q) . and PPMX dRoser et alJ |2008b cata- 



logues (Fig.|5]l. For UCAC2, there appears to be a periodic sys- 
tematic variation. However, when comparing each of them with 
Tycho-2, no similar effect was detected, indicating that this prob- 
lem affects mostly faint stars, whose first epoch data for deriving 
proper motions were neither AC2000 nor Cartes du Ciel plates. 
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Fig. 2. Comparison between some results of this work and those 
published in the D06 catalogue. The differences in the open clus- 
ter proper motions are in the sense of hereof minus D06 results. 



For UCAC3, we measure strong offsets in both proper mo- 
tion components (4.3 + 1.6 mas yr in c and -3.1 + 1.2 mas yr"' 
in 6), while relative to PPMX there is a small offset in the right 
ascension component (1 .2 + 0.75 mas yr"'), but no noticeable ef- 
fect in declination. Since all these four catalogues use common 
material and cannot be considered to be completely independent, 
no firm conclusions could be drawn from a simple analysis and 
a detailed study is required to address the origin of these effects. 

For clusters classified as "non-existent", "doubtful" or "not 
found" in the D07 catalogue, the chosen parametrisation failed 
to provide a good fit (our visual classification of the solutions 
was "poor"). This may indicate that those regions probably are 
not composed by two different kinematic components, thereby 
confirming the flag "non-existent" of the D07 catalogue. This is 
important, since the aforeme ntioned flag is based mainly on a vi- 
sual inspection performed bv lSulentic & Ti fft ( 1979). These ob- 
jects are the NGCs 2224, 2234, 2265, 6525, 673^1 6837, 6839, 
6840, 6843, 6858, 6950, and 7084. We note that this "poor" clas- 
sification cannot be caused by the limiting magnitude of our cat- 
alogue since NGC objects should be relatively bright. The clus- 
ter Dol Dzim 2 and Dolidze 35, with the flags "not found" and 
"doubtful" in the D07 catalogue were also classified as present- 
ing a "poor" solution. 

Regarding all the other clusters, we obtained solutions classi- 
fied as "good" for the NGCs 1817, 2169, 2194, 2355, 2682, and 
Berkeley 82. The clusters for which we obtained solutions con- 
sidered as "intermediate", and therefore require further analysis, 
likely by using additional information from photometric, radial 
velocity data and deeper or more precise proper motions, are: 
Berkeley 43, Berkeley 45, Berkeley 47, the Chupina series. King 
26, the NGCs 1807, 2678, 7036, Roslund 1, Skiffs J0614+129, 
and Teutsch 1 1 . 



' Confirming the results of lBoeche et all ( 120031) , which was based on 
Tycho-2 data. 
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Fig. 3. Systematic effect in the results for the common clusters 
between this work and the D06 catalogue. The dashed line is a 
weighted least-squares fit on the points. 



Fig. 4. Differences in proper motions between common PM2000 
and UCAC2 stars in the sense PM2000 minus UCAC2. This 
comparison includes only stars in the regions of the common 
open clusters reduced by this work and D06. The dashed line is 
a least squares fit to the points. 



In the case of NGC 2678, we note that the stars classified as 
members are concentrated in direction of NGC 2682. The bright- 
est stars in this region (those that visually define the cluster) were 
also not classified as members. Since these two clusters are close 
to each other, this may indicate that these stars are members of 
NGC 2682 and that NGC 2678 does not exist. 

The remaining objects on our list were classified as "poor". 
In the case of Berkeley 29, even though we have around 80% of 
the stars with membership probability greater than 5 1 % in com- 
mon between our and D06 analyses, the cluster's proper motions 
in right ascension are not compatible with those published in 
D06. We notice that both solutions are probably only good fit- 
tin gs of the field d istribution, since from the photometric work 
of Tosi et all (12004) we can see that the vast majority of this 
open cluster's stars are at V > 18, far beyond the reach of the 
PM2000 (used here) or UCAC2 (used in D06) catalogues. For 
NGC 2395, we could not distinguish a concentration of member 
stars, and the bright stars that define the cluster were not consid- 
ered as members. 

For all the other solutions classified as "poor", we exam- 
ine with regions in which the difference between the number of 
field and the cluster stars was very high. As a result, the nat- 
ural tendency of the likelihood function was to reach a better 
fitness value for the determination of a field only distribution. 
Nonetheless, we note that some clusters classified as "poor" can 
have a good determination of their proper motions, as in the case 
of NGC 7772, which we comment on the next section. 



5.1. Comments on selected objects 

Because of the physical interest of some open clusters in this 
study, we briefly discuss the properties of a small subset of them 
individually. 



5.1.1. NGC 1807 

The NGC 1807 open cluster is a concentration of stars close 
to NGC 1817. The solutions obtained for this object using all 
the four parametrisations indicate that its most probable mem- 
bers are concentrated in one border of the extracted field, mainly 
in the direction of NGC 1817. For the 4D parametrisation, the 
concentration in that part of the star chart is even larger than 
the other ones. In addition, the fitting parameters obtained for 
the cluster's PDF are compatible (within the fitting errors) with 
those obtained for NGC 1817 (as can be seen in Table HI . 

We conclude that the existence of a separate cluster 
from NGC 1817 in the NGC 1807 region is not supported 
by our analy sis of P M2000 kinematical data. This upholds 
Balaguer-Nunez et alJ ( [2004a.b) assertion, that there is actually 
only one extended cluster in that region. However, we note that 
a radial velocity study would be of great value to definitely settle 
this issue. 



5.1.2. NGC 2194 

NGC 2194 is an open cluster located in a moderately rich field 
projected along the Galactic anti-cen t re direction , whic h has 
been sfiidied bv iKveong et alJ (l2005h . iPiatti et alJ (l2003l) . and 
ISanner etal ](j2000^ Only the final of these three aforementioned 
studies deals with proper motions. Sanner et al. (2000) were the 
first to analyse this object's proper motion. In their study, proper 
motions were obtained by using photographic plates from the 
Bonn Doppelrefraktor (first plate from 14.02.1917) and from 
CCD observations performed at the 1 .23m telescope at the Calar 
Alto Observatory (15.10.1998). They determined the cluster's 
proper motion to be (//„ cos 6; fig) = (-2.3 + 1 .6; 0.2 + 1 .5) mas 



yr (the indicated errors are the fitted widths of the PDF). 
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a(h) a(h) 

Fig. 5. Average differences in bins of a in the proper motion components of common stars between PM2000 and UCAC2, UCAC3, 
and PPMX catalogues. 



This cluster was later analysed, automatically, in D06, in 
which its proper motion was computed from UCAC2 data 
(-0.59 + 0.53; -3.49 + 0.53) mas yr"'. As one can promptly 
see, there is a wide discrepancy between the proper motion val- 
ues obtained by those two studies, even if compatible to within 
3cr. It is even more confusing that the entry in the D07 catalogue, 
(-0.31 + 0.64; -4.40 + 0.64) mas yr"', is once again differenj^ 
even though it is compatible with the one obtained using UCAC2 
proper motions. 

We found herein the mean proper motion of this object to be 
more compatible with that of D06, with a value of (0.07 + 0.20; 
-2.69 + 0.29) mas yr ' . The almost 3cr discrepancy between the 
cluster's proper mo tion results obtained by this work, D06, and 
ISanner etal H2000), may be caused by the different material- 
izations of the HIPPARCOS reference system. As for UCAC2, 
PM2000 is l inked to HIPPARCO S as materialized by Tycho-2. 
In constrast.l Sanner et al.l (l2000l) used ACT stars to perform the 
transformation from plate to celestial coordinates, ACT being 
linked to HIPPARCOS, but as materialized by Tycho-1. 

It is also interesting that the proper motion of this cluster 
is compatible with that obtained for Skiff J0614+129. We can 
also identify some of the brightest stars in the region of Skiff 
J0614+129 (those that define the central concentration in its 
member star chart) among the list of NGC 2194 members, but 
a more detailed analysis using photometric and radial velocity 
data would be required to check whether they are somehow con- 
nected. 



^ The same value is quo ted in the WEBDA database, and was com- 
puted from Tycho-2 data in lDias et 



5.1.3. NGC 7772 

NGC 7772 was fi rstly proposed to be an open cl uster remnant by 
iBica et alj (1200 Ih . and confirmed to be one by ICarrarol ( |2002|) . 
In the latter study, NGC 7772's age was estimated to be 1.5 Gyr, 
and its distance from the Sun to be 1 .5 kpc, which assigned the 
membership to 14 stars. Nonetheless, since that work was purely 
based on photometric data, Carraro himself warned that his anal- 
ysis needed to be constrained by radial velocity and proper mo- 
tion determinations. Here we constrain the membership of the 
stars in this object by using proper motions. 

We note that the fitting of two kinematic populations to this 
cluster VPD, using the PDF that does not take into account the 
internal dispersion, does not allow a highly reliable segregation 
between the remnant and the field, since only two stars were 
classified as members using the criteria for the number of mem- 
bers presented at the beginning of this section (n = 2.34 + 1 .56). 
We believe that this to be caused by our choosing the PDF that 
does not take into account the internal dispersion between the 
member stars: since this is a remnant, the dispersion could be 
greater than the individual measurement errors, and we should 
be able to resolve it. 

Since this is a sparsely populated region, we can, however, 
analyse a plot of the proper motions represented as vectors su- 
perimposed directly on the Star Chart, as can be seen in Fig.|6] 
This graphical representation allows us to observe what is hap- 
pening from the point of view of the kinematic, spatial, and mag- 
nitude distributions. We realise that there exists a tightly packed 
spatial concentration of similar magnitude stars that share almost 
the same proper motion vector in this region. 

On the basis primarily of their proper motions, we visually 
chose the four (compatible with n to within ~ Icr) most simi- 
lar stars in the high-concentration region, two of which are the 
member stars classified by our automatic reduction. We note that 
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Fig. 6. Star chart of the field around NGC 7772. The non- 
member stars have their proper motions represented by the grey 
vectors, and the members by black ones. Meridian V magnitudes 
are also shown. 



the other stars have higher membership probabilities in the 4D 
solution. 

Since these are the constituent stars of this object, and whose 
errors in their individual proper motions are very small ((///e^) - 
4%), the average proper motion is a very good estimate of this 
remnant's proper motion. The final average proper motion of the 
object is therefore (/i„cos(5> = 14.90 + 3.13, (fig) - -10.65 + 
1.11 mas yr which is highly compatible with that determined 
automatically (see Table |4|. This is the first ever determination 
of a proper motion for this object. 

5.1.4. Berkeley 82 

Berkeley 82 is a sparsely populated open cluster that has yet to 
be thoroughly studi ed. This cluster is projected in the direction 
of a HI supershell (iKim & Kool 12000 ) that has a kinematically 
derived distance of 1 .4 kpc and a projected surface of 340 x 540 
pc^ (in the 1, b directions). These authors suggested that this ~ 5 
Myr supershell may be physically related to some of the open 
clusters in its region: NGC 6738, Berkeley 43, and Berkeley 82. 
The existence of the first of these three objects was ruled out. 
Berkeley 43 perhaps is an interesting object without proper mo- 
tion determination in the literature, but unfortunately our auto- 
matic solution for it was classified as "intermediate" probably 
because of the faintness of this cluster. 

However, for Berkeley 82, we were able to assign member- 
ship of stars to this object and compute its proper motion, which 
was measured to be (3.02 ± 0.61; -1.64 + 0.48) mas yr"'. From 
D07, a Tycho-2 determinati on for the pro per motion of this ob- 
ject was obtained by Lokti n & Beshenovl (2003) during a study 
of the rotation rate of the Galaxy (-0.06 + 0.99; -4.38 + 1.78) 
mas yr Nonetheless, no list has been published with individ- 
ual membership probabilities, so one cannot explore the connec- 
tion between the cluster and the shell on a star by star basis. 



According to WEBDA Berkeley 82 is located at a distance of 
870 pc, indicating that this cluster could well lie at one of the 
borders of the shell (if the cloud size is taken into consideration). 

Based on the quoted distance, our proper motion determi- 
nati on amounts to a t angential velocity of ~ 13.8 + 2.4 km 
s~'. lKim & Kool (l2000l) . determine the expansion velocity of the 
shell to be ~ 15 km s"', and its central velocity to be 18 km s"'. 
The similarities between the velocities that they obtained for the 
HI supershell and that which we obtained for the open cluster are 
quite remarkable. A radial velocity study would undoubtedly be 
of great value to determine the space velocity of this object, and 
finally reach a conclusion about the physical connection between 
this cluster and the supershell. 

5.1.5. NGC 1663 



During their photometric analysis, ' Baume et al.l (l2003h sug- 
gested that NGC 1663 is a ~2 Gyr possible open cluster remnant. 
In that study, the authors derived a list of 2 members and 4 possi- 
ble members, from a classification based on UBVI photometric 
data. They noted, however, that it was difficult to determine the 
true nature of this object. 

Our automatic solution for NGC 1663 was classified as 
"poor", since we could observe neither a clear concentration of 
member stars, nor any correlation with the brightest stars, in the 
member star chart. Moreover, the member list obtained by our 
solution included only one star from the photometric member 
Ust. 

Unfortu nately, it is unclear whether members and possible 
members in Baum e et al. I (l200l show common proper motions. 
Because of the lack of ancient epoch data in this particular zone 
of PM2000, the proper motion errors are greater than the average 
of the catalogue. This lack of precision prevents us from drawing 
firm conclusions, and a more precise proper-motion study and/or 
a radial velocity analysis is needed to clarify the nature of this 
possible open cluster remnant. 

6. Conclusions 

We have developed a fully automatic system to determine the 
kinematic parameters and membership probability lists of open 
clusters using kinematic and spatial data. The adopted method 
permits us to take into account any additional parameter (such 
as radial velocity), given that an analytical form of its probabil- 
ity distribution functions for the cluster and field populations are 
known a priori. Using this tool, we visually compared the re- 
sults for several possible parametrisations for the PDFs for all 
open clusters in PM2000. We concluded that the most reliable 
function is one that does not take into account the intrinsic dis- 
persion of the cluster. 

Based on PM2000 data, we obtained proper motions and 
kinematic membership lists for open clusters in the Bordeaux's 
PM2000 zone. For five of them, it was the first such measure- 
ment in the literature. We note, however, that for some of those 
objects, additional studies are necessary, if possible, based on 
radial velocity data to confirm the membership determination. 
Moreover, we confirmed from the kinematic point of view, the 
non-existence of 13 NGC objects, as well as Dol Dzim 2 and 
Dolidze 35, which are flaged as "not found" and "doubtful" in 
the D07 catalogue. 

We concluded that the open cluster NGC 1807 proba- 
bly does not exist, as previ ously suggested in the literature 
dBalaguer-Nunez et al.ll20043lbl) . Nonetheless, we note that a ra- 
dial velocity study to finally settle this issue is still required. 
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We also determined the tangential velocity of Berkeley 82 
to be remarkably similar to the velocities associated with the HI 
supershell located near it. Thus a physical connection between 
both objects remains plausible. 

By comparing PM2000 and UCAC2 proper motions, we 
found a periodic systematic variation as a function of the object's 
right ascension. Strong offsets in both proper motions compo- 
nents were also detected on comparing PM2000 with UCAC3. 
However, a comparison of PM2000 with PPMX showed only a 
small offset in the fia cos 6 component, and no noticeable effect 
in the fig component. 

We finally conclude that the blind use of parametric fully au- 
tomatic methods in present large-scale proper-motion catalogues 
is unreliable for most open clusters, since incorrect conclusions 
are often derived by inspecting VPD and probability histogram 
plots alone. When conducting studies based on these types of 
analyses, the results should be critically assessed using indepen- 
dent data; one should avoid using the aforementioned diagrams, 
since they are directly connected to the fitting itself. As a re- 
sult, one only checks how mathematically successful the fit is, 
and not its physical reality. The use of additional data, such as 
readily available star charts (to identify the spatial clustering of 
member stars), is highly advisable in the absence of multi-colour 
photometry or large-scale radial velocity data. 
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Table 1. Our cluster sample sorted by class type. The parameter values are taken from the D07 catalogue at Vizier (iDias et al 
apart from the coordinates of the clusters Berkeley 29, 43, 45, and 47 since they were re-centred by means of a detailed visual 
examination performed with Aladin (*). The Galactic coordinates presented here were computed by VizieR. Classes: A-reference, 
B-known, C-known without proper-motion and D-others. In the remarks, pocr indicates the possible open cluster remnants while the 
"c." indicates which condition of the section 2 .2 applies. Note: the Chupina clusters are actually subclusters formed on the outskirts 
of NGC 2682 (IChupina & Vereshchaginll998h . 
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Table 4. The results obtained for all the extracted fields in the PM2000 zone. The c subscripts indicate the parameters for the clusters, while the/ subscripts indicate the fields ones. The correlation 
coefficient (pf), mixture proportion (n^) and estimated number of members (iimemb) are adimensionals, whereas all the other parameters are in mas yr"' units. The "Class" column is the same as in 
Table[T] The solutions are divided following their visual classifications. 
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0.13 


78 
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0.21 


-3.59 


+ 


1.02 


-3.17 


+ 0.95 
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17 
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0.15 


0.23 


+ 


0.08 


12 


B 
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+ 


0.35 
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+ 


0.31 
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+ 
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5.05 


+ 


0.81 
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+ 


0.87 


0.07 


+ 


0.19 


0.32 


+ 


0.11 


15 


B 


NGC 7036 


-5.04 


+ 


0.41 


-6.84 


+ 


0.35 


0.06 


+ 


1.26 


-8.15 


+ 1.11 


8.23 


+ 


0.98 
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+ 


0.93 


0.56 


+ 


0.09 


0.26 


+ 


0.08 


20 


B 
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7.18 
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0.89 
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+ 
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+ 


0.12 


0.47 


+ 


0.08 


70 


B 
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0.30 


+ 


0.78 
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+ 


0.68 


-0.36 


+ 


1.02 


-3.52 


+ 1.65 


4.38 


+ 


1.02 


8.13 


+ 


1.31 


0.08 


+ 


0.16 


0.36 


+ 


0.15 


23 


C 


Berkeley 45 
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+ 


0.96 
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+ 


0.94 


6.28 


+ 


1.13 


-3.94 
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+ 


1.17 


3.75 


+ 


1.10 


-0.37 


+ 


0.33 


0.47 


+ 
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42 


C 
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0.90 
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0.86 
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+ 
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+ 2.46 


6.79 


+ 


1.63 


9.65 


+ 


2.05 


0.65 


+ 


0.14 


0.60 


+ 


0.12 


35 


C 


King 26 


3.80 


+ 


0.58 


-5.72 


+ 


0.47 


1.18 


+ 


1.99 


-6.31 


+ 1.73 


7.68 


+ 


1.72 


6.50 


+ 


1.58 


0.04 


+ 


0.22 


0.45 


+ 


0.15 


21 


C 


Skiff J0614+129 


-0.31 


+ 


0.29 


-2.73 


+ 


0.37 


-0.57 


+ 


0.44 


-4.31 


+ 0.63 


4.48 


+ 


0.32 


6.56 


+ 


0.47 


-0.05 


+ 


0.09 


0.18 


+ 


0.05 


29 


C 


NGC 1807 


2.06 


+ 


0.23 


-4.56 


+ 


0.24 


3.57 


+ 


0.68 


-7.83 
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7.80 


+ 


0.57 


8.35 


+ 


0.59 
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+ 


0.07 


0.66 


+ 


0.03 


529 


D 


NGC 2678 


-1.94 


+ 


0.32 


-2.00 


+ 


0.28 


-4.18 


+ 


0.75 


-5.34 


+ 0.57 


8.29 
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0.55 


6.10 


+ 


0.43 
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+ 


0.09 


0.11 


+ 


0.04 


16 


D 
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-2.24 
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0.30 
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1.28 
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16 


D 


Poor 


Berkeley 29 
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0.55 
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Table 4. continued. 
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Table 5. Example of the membership probability table for the remnant NGC 7772. The PM2000 and 2MASS columns points to the identifiers of the object in the respective 
catalogues. The Prob. column indicates the membership probability as calculated from eq.|9] 
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